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Résumé 

La stabilité d'une variable aléatoire peut être généralisée dans un 
cône convexe. Les résultats concernant la représentation de LePage et 
le domaine d'attraction sont analogues mais différents de ceux qui sont 
obtenus dans l'espaces de Banach. Quelques exemples de lois stricte- 
ment max-stables et a-stables sont présentés dans cet article afin de 
montrer le lien entre le processus ponctuel poissonnien et la loi stable 
dans un cône convexe. 

Mots clés : a-stable, max-stable, cône convexe, processus ponctuel poisson- 
nien, série de LePage, simulation. 



1 Introduction 

Un vecteur aléatoire X à valeurs dans un espace de Banach a la loi stricte- 
ment a-stable {StaS) si pour tous a,b > 

a^/"Xi + 6i/"X2 = (a + 6)i/"X, (Ll) 

OÙ Xi,X2 sont les copies indépendantes de X, et = représente l'égalité en loi. 



La stabilité ( 1.1 ) peut être définie dans tous les espaces où l'addition des élé- 
ments et la multiplication par les scalaires positifs sont bien définies, i.e. dans 
un cône convexe. Dans ce cas les résultats principaux (la représentation de 
LePage, le domaine d'attraction, etc.) sont analogues mais différents de ceux 
qui sont obtenus dans l'espace de Banach. L'intervalle des valeurs possibles de 
l'indice caractéristique a dépend des propriétés algébriques et topologiques 
de base de l'espace où l'élément aléatoire StaS est défini. Par exemple il est 
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connu que a est dans l'intervalle (0, 2] pour les lois a-stables dans l'espace de 



Banach. Si on remplace l'addition dans (1.1 ) par les opérations de maximum 
et minimum, les vecteurs aléatoires max-stables et min-stables apparaissent. 
Les lois max-stables existent pour tout a > 0, tandis que les lois min-stables 
peuvent avoir toutes les valeurs négatives pour l'indice a. Davydov et ses co- 
auteurs ont caractérisé dans |3| les lois des éléments aléatoires StaSs dans 
un cône convexe et présenté les valeurs possibles pour l'indice a selon les 
propriétés du semigroupe et la distributivité de l'opération de multiplication 
correspondante. Nous présentons dans cet article quelques exemples des lois 
StaSs et max-stables dans avec la simulation numérique de la densité, 
afin d'illustrer les rôles de l'indice a et de la mesure spectrale dans la forme 
de la densité de loi StaS dans différents cônes convexes. 

Rappelons qu'un vecteur aléatoire a-stable symétrique ou avec a < 1 à 
valeurs dans admet la représentation de LePage suivante [9] 



fc=l 



(1.2) 



où c > 0, Ffc représente la somme de k variables i.i.d. de loi exponentielle 
d'espérance 1 et ei, e2, • • • sont les vecteurs aléatoires unitaires i.i.d. et indé- 
pendants de {r^}^^. Notons IIq, le processus ponctuel dont le support est 

{cT~^^°'ek, k > 1}, i.e. 

CX) 

k=l 

OÙ ôx désigne la mesure de Dirac au point x. Le support du processus IIq 
est un ensemble fermé aléatoire union-stable \12\ ^4], i.e. pour tous a,b > 



a 



l/a 



supp(n;) U 6i/"supp(n'^) = (a + 6)i/°supp(n„). (1.4) 



où et H'^ sont les copies indépendantes de IIq. Cela signifie que l'union- 
stabilité du support de processus IIq est étroitement liée à la stabilité clas- 



sique (1.1). Cette relation est l'idée clé dans |3j pour obtenir les résul- 
tats concernant la représentation de LePage et le domaine d'attraction de 
loi StaS dans un cône convexe. Les exemples présentés dans cet article 
montrent que les processus ponctuels sont utiles non seulement pour obtenir 
les résultats généraux sur les lois StaSs dans un cône convexe, mais aussi 
pour déduire la fonction de répartition et la régularité de la queue d'une loi 
max-stable avec la mesure spectrale donnée. 

Un autre but de cet article est de montrer que les processus ponctuels 
poissonniens peuvent être utilisés pour simuler les vecteurs aléatoires suivant 
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les lois StaSs et max-stables. Ces deux familles de lois apparaissent souvent 
dans les applications du domaine des événements rares qui demandent le 
traitement de valeurs extrêmes. Il existe des méthodes connues pour engen- 
drer les données de loi stable |2] |23| |10| et de loi max-stable |,13j [22\ [21j . 
En utilisant la représentation de LePage généralisée dans un cône convexe, 
nous pouvons engendrer les données suivant les lois StaSs et max-stables de 
façon analogue. Notre méthode ne nécessite pas la discrétisation de la me- 
sure spectrale, elle est donc plus pratique pour simuler les vecteurs aléatoires 
StaSs et max-stables avec la mesure spectrale absolument continue. 

Nous résumons dans la section |2] les définitions et les résultats de base 
concernant les lois StaSs dans un cône convexe. Pour la présentation plus 
détaillée nous renvoyons le lecteur à l'article La section [3] présente les 
lois strictement max-stables et a-stables avec des exemples concrets. Nous 
obtenons l'expression explicite des lois jointes et marginales des vecteurs 
aléatoires strictement max-stables avec les mesures spectrales données. Les 
résultats de simulation des vecteurs aléatoires dans les deux familles sont 
présentés. 



2 Lois strictement stables dans un cône convexe 
2.1 Cônes convexes 

Définition 2.1. Un cône convexe IK est un semigroupe abélien topologique, 
supposé complet et séparable, avec une opération de multiplication par des 
scalaires positifs, (x, a) — )■ ax, continue pour tout x £ IK et a > tel que les 
conditions suivantes sont remplies : 

1 ) a{x + y) = ax + ay, a > 0, x,y G IK 

2) a{hx) = {ab)x, a,b> 0, x e IK 

3) lx = x, xe IK 

4) ae = e, a > 0, e est l'élément neutre de IK. 

Le cône IK s'appelle cône pointé s'il y a un élément unique dit l'origine 
tel que ax — )• lorsque o ^ pour tout x £ IK \{e}. 

Remarquons que la condition suivante qui est typiquement imposée dans 
la littérature (voir par exemple |i8.j) pour un cône n'est pas nécessaire ici 

(a + b)x = ax + bx, a,b> 0, x e K. (2.1) 

Par exemple M+ = [0, oo) avec l'opération d'addition définie par le maximum 



ne vérifie pas cette condition. Dans un cône pointé, si la distributivité (2.1) 
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a lieu ou s'il existe l'élément x ^ e qui possède l'inverse (— x), alors e = 0. 
Généralement l'élément neutre ne coïncide pas nécessairement avec l'origine. 



Définition 2.2. Un élément z & IK est a-stable avec a ^ 0, si 

= (a + 6)^/"z 

pour tous a,b > 0. 

Dans la suite on note IK(a) l'ensemble des éléments a-stables de K. Il 
est clair que e, G IK(a) pour tout a ^ 0. 

Définition 2.3. Le cône IK s'appelle un cône normé si IK est métrisable par 
une distance d qui est homogène à l'origine, i.e. d{ax,0) = a(i(x, 0) pour 
tout a > et X Çz IK. La valeur \\x\\ = d{x, 0) s'appelle la norme de x. 

Dans la suite on suppose que IK est un cône normé. Il est clair que ||x|| = 
si et seulement si x = 0. De plus si e 7^ 0, alors la propriété 4) de la définition 



2.1 implique que ||e|| = d(e, 0) = 00. Il est donc essentiel d'admettre que d 
peut prendre la valeur infinie. Par exemple, si IK est le cône M-|- = [0, 00] 
où l'opération d'addition est définie par le minimum, noté (M+, A), alors la 
distance euclidienne d'un élément non-zéro x G M+ à 00 (étant l'élément 
neutre) est infinie. 
L'ensemble 

S = {x \ \\x\\ = 1} 

s'appelle la sphère unité. Le cône pointé admet la décomposition polaire réa- 
lisée par la bijection x O (||x||, x/||x||) entre 

K' = IK\{0,e} 

et (0,00) X S. 

La distance ou la norme dans IK est dite sous-invariante si 

d{x + h,x) < d{h,0) = \\h\\, x,h £K. 

Cette contrainte n'est pas triviale, par exemple le cône (M-|-, A) avec la dis- 
tance euclidienne ne vérifie pas cette propriété. Si IK a une norme sous- 
invariante, alors = e, et pour tout a G (0,1), e est l'élément unique de 
norme finie qui appartient à IK(a). 

Les exemples typiques de cônes qui possèdent ces propriétés sont les 
espaces de Banach ou les cônes convexes dans l'espace de Banach, la famille 
des sous-ensembles compacts (ou compacts convexes) d'un espace de Banach 
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avec l'addition de Minkowski |5] [7] [6], la famille des ensembles compacts 
dans M'^ avec l'opération d'addition définie par l'union et la famille des 
mesures finies avec l'addition conventionnelle et la multiplication scalaire [3] 
|17| . Un autre exemple typique est l'espace = [0, oo) avec l'opération 
d'addition définie par le maximum, i.e. x + y = xV y = max(x, y). D'autres 
exemples et informations peuvent être trouvés dans |4]. 

2.2 Processus ponctuels dans un cône 

Considérons un cône normé IK muni de la u-algèbre borélienne B{¥^). La 
boule ouverte de rayon r centrée au point est notée par 

Br = {x£'K \ \\x\\ < r}. 

L'intérieur de son complément est noté par 

B'^ = {x \ \\x\\ > r}. 

Une mesure de comptage m est une mesure de forme suivante 

m = + -\ — = y^^x, 

i 

OÙ {xi,X2, ■ ■ •} est un ensemble de points dénombrable au maximum. Soit 
Mo (respectivement Ai) la famille des mesures de comptage m sur B(JK) 
telles que m{B^) < oo (respectivement m(Br) < oo) pour tout r > 0. 
Pour couvrir ces deux cas avec la même notation, on note au lieu de 
B^ (respectivement Br) dans le cas où l'on considère les mesures dans A^o 
(respectivement M). Alors nous avons toujours m{Ar) < oo, c'est-à-dire 
seulement un nombre fini de points Xi du support de m se situent dans Ar- 

Un processus ponctuel /i est une application mesurable d'un espace de 
probabilité à A^o (ou Ai) muni de la a- algèbre engendrée par les ensembles 
des mesures m £ Aio (ou m G Ai) telle que m{B) = n pour les ensembles 
boréliens C K et n > 0. 

Soit A une mesure sur IK qui est finie sur Ar pour tout r > 0. Un 
processus ponctuel H est un processus poissonnien de mesure d'intensité A, 
si pour toute famille d'ensembles boréliens disjoints Fi, . . . , Fn, les variables 
aléatoires Il{Fi), . . . ,Il{Fn) sont mutuellement indépendantes et distribuées 
suivant les lois de Poisson d'espérance A{Fi), . . . ,A{Fn) respectivement. 

Définition 2.4. (Processus ponctuel poissonnien stable Ila,a) Soit Aa^a une 
mesure sur IK' définie par rria x a où a est une mesure finie sur B{S) et nia 
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est une mesure sur (0, oo) telle que pour a ^ on a 



maiir, oo)) = r " si a > 0, ^ ^ Q 
ma((0, r)) = r""^ si a < 0, 

Le processus ponctuel poissonnien stable est un processus poissonnien de 
mesure d'intensité Aa,cr, notéIla,cr- 

La mesure a s'appelle la mesure spectrale. L'importance du processus 
vient du fait que le support de Ila,a est union-stable, i.e. Ila,a vérifie 
l'égalité ([L4|. 

Pour tout a 7^ 0, le processus Ila,a admet la représentation (1.3) (Th. 
3.3 |3j) avec les éléments aléatoires {e^} i.i.d. à valeurs dans S et de loi 
= a{-)/a{S) et c = C7(5)V". 

Soit {.Çfc, A; > 1} une suite d'éléments aléatoires i.i.d. à valeurs dans K. 
Un processus ponctuel empirique est défini par 



k=l 

OÙ {bn, n > 1} est une suite de constantes de normalisation telle que 

que pour avoir 



oo, si n — )• oo. On va montrer dans le théorème 2.5 



de propriétés intéressantes, les constantes 6„ doivent avoir la forme suivante 

6„ = n^/"L(n), n > 1, (2.3) 

avec a 7^ et L{n) une fonction à variation lente à l'infini. 

La convergence faible des processus ponctuels =^ A* a heu, si E/i(/i„) — )• 
E/i(//) lorsque n — )• oo pour toute fonction h sur M ou A^o à valeurs dans 
M continue dans la topologie vague et bornée. Le processus poissonnien Ila,a 
est la limite faible du processus empirique /3„, si la loi des a la queue à 
variation régulière [TB]. On donne ici ce résultat sous la forme de 

Théorème 2.5. (Th. 4-3 ^) Soient ^,^i,^2^ ■ ■ ■ des éléments aléatoires i.i.d. 
à valeurs dans IK' . Alors Pn =^ ^a,a lorsque n — )• oo pour a > si et 
seulement s'il existe une mesure finie a sur B{S) telle que on a 

lim nP G ||^|| > r6„ 1 = (T(5)r^", (2.4) 



n—^oo 



m 



pour tout r > et B G B{S) avec a{dB) = 0, où bn est de forme (2.3). 
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Remarque 1. La condition (2.4) est typique pour les théorèmes limites 



de somme des éléments aléatoires (voir |18) page 154 et [I] page 167). Il existe 
un résultat similaire pour a < (Cor. 4.4 

La loi d'un élément aléatoire ^ à valeurs dans IK ' est dite à queue régulière 



si la condition de variation régulière (2.4) est vérifiée. 



2.3 Série de LePage et domaines d'attraction de loi StaS 
Un élément aléatoire X dans IK a la loi strictement astable {StaS) s'il 



vérifie la condition (1.1) de l'introduction avec les opérations d'addition et 
multiplication définies sur K. Le théorème ci-dessous donne une famille riche 
de lois StaSs par leur décomposition en série. 

Théorème 2.6. (Th. 3.6 [40 Soient {Xk, k > 1} et {e^, A; > 1} deux suites 
indépendantes de variables aléatoires i.i.d.. Les variables aléatoires Afe ont la 
loi exponentielle standard, les éléments aléatoires sont de même loi à(-) qui 
est une mesure finie normalisée sur S. Notons = Ai-|-- • ■ + \k-, k > 1. Si la 
valeur principale de l'intégrale J xlla,à{dx) est finie, i.e. xlla,àidx) < oo, 
avec probabilité 1, alors pour tout z E IK{a) et c>0, la série 



^ + (2.5) 
fc=i 



converge presque sûrement (p. s.) et ^^,0- admet une loi StaS sur IK. Le para- 
mètre a s'appelle l'indice caractéristique et la mesure a{-) = c°(t(-) s'appelle 
la mesure spectrale. 



Si la norme de IK est sous-invariante, alors la somme infinie (2.5) converge 
absolument p. s. pour tout a G (0, 1). 

Définition 2.7. Un élément aléatoire X à valeurs dans IK appartient au do- 
maine d'attraction d'un élément aléatoire StaS, noté ^a,a > si pour {Xn, n > 
1}, une suite de copies indépendantes de X, on a 

b-\Xi + --- + Xn)^Ca,a, (2.6) 

OÙ {bn, n > 1} est une suite de constantes de normalisation positives et 
désigne la convergence faible des éléments aléatoires à valeurs dans IK. 

Le fait que l'élément aléatoire X appartienne au domaine d'attraction de 
(,a,a sera noté par l'écriture "X G Dom(Q;, a)". Le théorème suivant montre 
qu'avec une condition supplémentaire sur IK la régularité de la queue de la 
loi de X implique que X G Dom(a, a) pour < a < 1. 



7 



Théorème 2.8. (Th. J^.l Supposons que IK a la norme sous-invariante. 
Si la loi d'un élément aléatoire X £ IK' vérifie la condition (2.4) avec a G 
(0, 1), alors X appartient au domaine d'attraction d'un élément aléatoire 



StaS qui admet la représentation du type de LePage (1.2) avec Tj^, et c 



qui sont définis comme dans le théorème 2.6 



Remarque 2. La condition (2.4 ) est nécessaire et suffisante pour qu'un 
vecteur aléatoire à valeurs dans 



appartienne au domaine d'attraction 



d'une loi StaS avec < a < 2 |T]. 



L'existence de la mesure spectrale d'une loi StaS dans un espace de 
Banach séparable est connue. L'existence de la représentation de LePage 
d'un élément aléatoire StaS dans un semigroupe plus général est démontrée 
dans [4J. Dans la section suivante nous présentons quelques exemples de lois 
max-stables et a-stables. 



3 Quelques exemples et simulation 
3.1 Lois strictement max-stables 

On considère les lois stables dans le cône (M!^, V). Les vecteurs dans Mî|_ = 
[0,00)'^ sont notés par x = [xi, . . . ,Xd)- Les relations et opérations sont 
définies par composante, c'est-à-dire pour x, y G 

x < y signifie Xi < yi,i = 1, . . . ,d, 

X ^ y signifie Xi ^ yi,i = 1, . . . ,d, 

x + y := X Vy = (xi Vyi,. . . ,3;^ Vyd), 

ax := {axi, . . . ,aX(i), a > 0. 

Les rectangles sont notés par 

(a, b) = {(xi, . . . ,Xd) \ ai < xi < bi, . . . ,ad < Xd < bd}, 

[a, b] = {(xi, . . . ,Xd) \ ai ^ xi ^ bi, . . . ,ad ^ Xd ^ bd}. 

Prenons la Loo-norme, i.e. ||x|| = max(xi, . . . , Xd)- La sphère unité dans cette 
norme {x | ||x|| = 1} est notée par S!^~^. 

Un vecteur aléatoire X = (X^^) , . . . , X^*^^) a la loi strictement max-stable 
si pour tous a,b > 
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où Xi, X2 sont les copies indépendantes de X. Soit ei,e2, • • • une suite de 
vecteurs aléatoires i.i.d. à valeurs dans S^"^ suivant la loi a. Soit Ai, A2, . . . 
une suite de variables aléatoires i.i.d. de loi exponentielle standard et Fj = 
Yl]=i ^j- Les suites {ej} et {Fj} sont indépendantes. La série suivante converge 
absolument p. s. pour tout a > 



\Jr;'/''e,. (3.1) 

i=l 



Le résultat direct du théorème 2.6 indique que cette série produit un vecteur 
aléatoire strictement max-stable avec l'indice a > et la mesure spectrale 
a. Dans la suite on utilise l'écriture "A^5(^(a, a)" pour indiquer la loi stric- 
tement max-stable d-dimensionnelle d'indice a et de mesure spectrale a. 

Soit Ç Q,^CT un vecteur aléatoire de loi A^5rf(a,fT). Puisque ^a,a et la série 
(3.1 ) ont la même loi, le maximum partiel suivant donne une approximation 
de la loi de £,a,a 

i=l 

Notons 

€i = {ef\. . . , ef et tj = mm{i \ e^^ = 1}, j = l, . . . ,d, i = 1,2, . . . . 
Supposons que 

P{eS^') = 1} = > 0, j = l,...,d. 

Le vecteur aléatoire converge vers (,a,a en probabilité lorsque k tend vers 
l'infini. En effet 

d d 
+ io^A = P{3j tel que r, > k} <Y, P{rj > k} = ^(1 - pj)\ 

La dernière égalité vient du fait que la loi de tj est géométrique 

P{r, =n} = P{e? < 1, . . . , e^^li < 1, e^f ) = 1} = (1 - 
Notons le processus ponctuel 

00 

na,.(-) = ^5r-/.^C)' (3.2) 

i=l 
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où Fj et €i sont définis comme dans (3.1 1. D'après la représentation du pro- 
cessus ponctuel poissonnien stable (1.3) et la définition (2.2), la mesure d'in- 
tensité de Ila,a est rrioi x a. Remarquons que pour tout r > 0, puisque 
ma{{r, oo)) = < oo, on a 

P{n„,^({x| ||x|| > r}) < oo} = 1. 

Il y a seulement un nombre fini de points du support de Ila,a qui sont grands. 
Le fait que Çq^q- et la série (3.1) ont la même loi implique que la loi de Çq^o- 
ne dépend que des points les plus grands du processus IIq, o-. La loi de ^a^a 
est facile à calculer, pour x > 

PUa,a e [0,x]} = P{n„,,([0,x]'^) = 0} = exp(-m« X ct([0, x]^^)). (3.3) 



La représentation de LePage (3.1 ) peut être écrite sous la forme suivante 



V ^fe ' • • • ' V k j — y'-Ti ' • • • ' 



l/a(d) 



l/a^ 



La loi marginale de ^a,a peut être obtenue en utilisant le processus ponctuel 



j = 1, 



(3.4) 



1=1 



Si la loi de la composante e["'\ notée u'-'-', est connue, alors 



(i) 



PUS < Xj} = P{n,,,(.)((x,,oo)) = 0} = exp(-m, x ^(^^((x,-, oo))). 



(3.5) 



Dans les trois exemples ci dessous nous utilisons le lien entre le processus 
ponctuel poissonnien et la loi strictement max-stable pour déduire la fonction 
de répartition de la loi avec la mesure spectrale donnée. L 'exemple [T] présente 
la mesure spectrale discrète et concentrée sur les points d'intersection des 
axes et de la sphère unité. Dans ce cas le vecteur aléatoire strictement max- 
stable a ses composantes indépendantes. L'exemple [2] présente la mesure 
spectrale uniforme. L'expression explicite de la loi jointe est obtenue pour 
d = 2. L'exemple [3] présente la mesure spectrale dans concentrée sur deux 
points qui ne sont pas nécessairement sur les axes. L'exemple [T] pour d = 2 
est un cas particulier de l'exemple [3] Des vecteurs aléatoires strictement 
max-stables à valeurs dans avec les mesures spectrales présentées sont 
simulés en utilisant le maximum partiel de la représentation de LePage ^fc. 
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Exemple 1. Considérons le vecteur aléatoire ^a^a de loi M.Sd[oi^a) 



avec 



i=l 

Où e, = {e^\...,ef^), = ' = ^ et f = 1, > 0, i = 

t ' î 1^ Sinon 

1, . . . C'est-à-dire la mesure spectrale a concentrée sur les points où les 
axes et la sphère unité se croisent. 

Pour chaque i, e\ est la variable aléatoire de loi de Bernoulli de para- 
mètre Pi. En appliquant (3.5), on a la loi marginale de ^OL^a 

P{Ç« <x,} = exp(-Kx-"). 



En appliquant (3.3), on a la fonction de répartition de Ça^o-, pour x > 



PUa,a e [0,x]} = exp ^- YJ'^^^'^ 



En conclusion, le vecteur aléatoire strictement max-stable ^a,a a ses com- 
posantes indépendantes si la mesure spectrale est concentrée sur l'intersec- 
tion des axes et de la sphère unité. C'est un résultat connu qu'on peut trouver 
par exemple dans [18] (Cor. 5.25). La figure [ï] (a.) présente les courbes de 
niveau de la densité de loi strictement max-stable bivariée d'indice a = 0.75 
et mesure spectrale a{-) = 0.55|(i o)}(-) + 0.55|(o,i)}(-)- 

Exemple 2. Soit a la mesure uniforme sur 5^"^. Soit e le vecteur aléa- 
toire unitaire de loi a. Avant de considérer la loi A^5rf(a,(T), nous montrons 
que la loi de composante e*^*), i = 1, . . . , d, est 

^^'^{■) = ^I{i}(-) + '^/^(•)l6([o,i))(-), (3.6) 
où n est la mesure de Lebesgue sur [0, 1). Puisque ||e|| = 1, on a 

P{3i G {1, . . . , 4 tel que e^^^ = 1} = 1. (3.7) 
L'absolue continuité de la loi de e implique que 

P{e» = eO) = j} = 0. (3.8) 



En considérant (3.7), (3.8) et l'uniformité de la loi de e, on a 

P{e« = 1} = ^, i = l,...,d. (3.9) 
a 
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Ainsi on a 



P{e»G[0,l)} = '^, i 



,d. 



(3.10) 



Les égalités (3.9 ), ( |3.10 ) et l'absolue continuité de la loi de e impliquent ( |3.6| ). 

Soit (^a,a le vecteur aléatoire de loi MSd{o:,a) avec a uniforme. En ap- 
pliquant (3.5), on a la loi marginale de ^a^a- 

P{dl < X,} = exp 



a + d 



{a + l)d ' 

En effet en notant r/ la variable aléatoire de loi ma, on a 



ma X cjW((xi,oo)) = P{7/e^*^ > Xi} 

.(0 



P{7? > Xi, e):> = 1} + P{r?e^*^ > x^, e'^> £ [0, 1)} 



a + d 
(a + l)d^^ 



ai' 



d-1 



-dtdx 



De manière analogue on peut déduire de (3.3) la fonction de répartition de 
(,a,a- Il n'y a pas d'expression générale pour tout d. Dans le cas où fi = 2 la 
loi jointe de Çq^o- est 

F^^.Àxi,X2) = exp (^- ^(^Jj^ 1) ("^(1)^(2°"^ + (a + 2)X("^p^ , 

où X(x) = xi /\x2 et X(2) = 3;i V^2- La figure [l| (b.) présente les courbes de 
niveau de la densité de loi strictement max-stable bivariée d'indice a = 0.75 
et de mesure spectrale uniforme. 

On discute finalement la régularité de la queue de cette loi. En prenant 
l'ensemble B = [(1,0), (l,t)] on B = [(0, 1), {t, 1)],0 ^ t ^ 1, on a pour tout 
r > 

h = P{n„,,((r, oo) X B) / 0, UaAitr, oo) x B^) = 0} 

< PU„,, G (r,oo) X S} 

< P{n„,,((r,oo)xS)/0} 
= h, 



ou 

h -- 



;i-exp(-CT(B)r-"))exp(-a(#)(tr)-") et h = l-exp{-a{B)r- 
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On en déduit 

P{^a,cT G X (r, oo)} ~ a{B)r^'^ lorsque r — )■ oo. 



Donc la loi de o- vérifie la condition de variation régulière (2.4). 



a. " 







0.2 0.3 0.4 



0.6 0.7 O.S 0.9 



Figure 1 - Courbes de niveau des densités réelles (gauche) et empiriques 
(droite) de la loi A^52(0.75, o"), où a est discrète définie dans l'exemple [ï] 
(a.), uniforme définie dans l'exemple |2] (b.). Le nombre de données simulées 
est 10^. 

Exemple 3. Considérons le vecteur aléatoire S,a,a de loi MS2{ot,(T) 
avec a concentrée sur deux points 



(3.11) 



OÙ ei = (1, a), 62 = (6, 1), < a, 6 < 1 et p + g = 1. Il est facile de voir que 
è,a,a est à valeurs dans la région entre deux demi-droites sortant du point 
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d'origine vers les points ei et 62- On en déduit la fonction de répartition 



exp(— px-^ " — a < X2/X1 < 1/b, 

P{Ç G [0,x]} = <( exp(-(pa" + g)x2"), X2/X1 < a, 
exp(— + X2/X1 > 1/6. 

En notant le moment d'arrêt r = min{z | ej = 62}, considérons la probabilité 
que S^a,a est à valeurs dans la demi-droite {ei x (0, 00)} 

00 

P{Ça,<xeeix(0,r)} = 5^P{ei=ei,r = i,r-^/"<ar-^/",r-^/"<r} 

1=2 

i-1 TD / ^2 H h Aj 

^ '^'^1 

1=2 ^ ^ 



> a " - 1, Al > r" 



pq{l - a")\ / / , -as -a^ 

P ^ ^ exp(-(p + ga )r ) 



De manière analogue on a 



P{Ça,.Ge2X (0,r)}= (^g-^^^-^Jexp(-(g + p6->-°). 

Si a = 6 = 1, alors p = g = 1 et la loi de S^a,a est définie sur la demi droite qui 
divise le premier quadrant en deux moitiés égales. Si a = 6 = alors S,a,a a 
ses valeurs strictement à l'intérieur du premier quadrant et ses composantes 
sont indépendantes. Cela coïncide avec l'exemple [î] Plus a et 6 sont proches 
de 1, plus il est possible que ^q^o- est à valeurs sur les demi-droites. La figure|2] 
présente la dépendance des paramètres a, 6 et a de la densité de AiS2{o(, cr) 
avec a définie par (3.11). 

Considérons la régularité de la queue de cette loi. Définissons deux vec- 
teurs aléatoires 



^62- 



Notons e\ , €2 \ ■ ■ ■ les copies indépendantes de e^*', i = 1,2. D'après le ré- 
sultat de l'amincissement du processus (voir par exemple Prop. 4.4.1 |19)). 
les processus IIi et 112 définis par 

00 

Ili = (5 -i/c (i) , i = l,2, 

. -, i j 

sont indépendants et poissonniens. Les mesures d'intensité sont respective- 
ment nia X pSeï et ma X qôe2- Prenons l'ensemble B = [(1, s), (1, t)], < 
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s<a<t<l. Il est facile de voir que pour tout r > 

Il = P{ni((r, oo) X {ei}) / 0, n2{{tr, oo) x {ez}) = 0} 

< PU„,, G (r,oo) x5} 

< P{^i((r,oo) X {ei}) /O} 
= h, 



ou 



/i = (1 — exp(— pr"")) ex.p{—q{tr)~°') et /2 = 1 — exp(— pr~"). 

Les deux termes Ii et I2 sont équivalents à pr~°' lorsque r — s- 00, ainsi on a 

P{Ço,CT & B X (r, 00)} ~ pr^" lorsque r — )• c«. 

De même manière, pour B = [(s, 1), (t, 1)], < s < 6 < t < 1 on a P{^q,^o- G 
(r, 00) X i?} ~ qr~°' si r — )• 00. La loi de Çq^o- vérifie la condition de variation 



régulière (2.4). 



3.2 Lois a-stables 

La difficulté dans l'étude de lois a-stables est que, sauf pour les cas parti- 
culiers (a = 0.5, let2), il n'y a pas de forme explicite pour la densité. Les 
vecteurs aléatoires a-stables sont déterminés par leur fonction caractéris- 
tique. 

Définition 3.1. La fonction caractéristique d'un vecteur aléatoire a-stable 
(0 < a < 2) dans s'exprime par l'expression suivante : 

0a,a(i) = exp ( --^ / 'il;a{{t,s))a{ds) + i{ô,t] 

où S'^^^ = {x G M"' I ||x|| = 1}, a est une mesure finie sur la sphère unité 
S'^~^, 6 est un vecteur dans Mf^, {t, s) représente le produit scalaire, 



l-a 



a / 1, 



C^ = \ r(2-a)cos{W2)' (3.12) 

l ï> a = l, 

et 

V'q(x) 



\x\ ( 1 ~ ^ sign{x) tan — j , a 1 



1^1(1 + sign{x) In |x| ) , a 
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Figure 2 - Courbes de niveau des densités réelles (gauche) et empiriques 
(droite) de la loi A4S2{a,(j) où a est discrète définie dans l'exemple |3] avec 
p = q = 0.5, Q = 0.75 et a = 6 = 0.2 (a.), a = 0.75 et a = 6 = 0.5 (b.), 
a = 3eta = b = 0.5. Le nombre de données simulées est 10^. 
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La loi stable dans M'^ est spécifiée par un nombre a entre et 2 dit indice 
de stabilité, une mesure finie a sur S'^~^ dite mesure spectrale et un vecteur 
6 dans M"^. Il existe des méthodes numériques pour obtenir la densité et la 
fonction de répartition approximatives d'une loi stable [15j . 

Dans le cas unidimensionnel la sphère unité ne contient que deux points, 
i.e. = { — 1, 1}. La mesure spectrale a se réduit à deux valeurs cr{ — l) et 
cr(l). La fonction caractéristique peut être écrite comme suit 




(^1 - î/3sign(t)tan^^ + iôt^ , a / 1, 
--^\t\ (l + sign(t) In \t\j + iôt^ , a = 1, 



où 7 = cj(l) + cr( — 1) et (3 = {a{l)—a{ — l))/j. La loi stable unidimensionnelle 
est déterminée par quatre paramètres : a G (0, 2], /3 G [—1, 1], 7 > et G M. 
Le paramètre /? contrôle l'asymétrie. La loi est dite biaisée totalement vers 
la droite si /3 = 1 et biaisée totalement vers la gauche si/3 = — 1. Si/3 = 
la densité est symétrique par rapport à ô. Les paramètres 7 et (5 sont les 
paramètres d'échelle et de position. Nous utilisons la notation S(i{a,a,ô) 
pour la loi a-stable d-dimensionnelle avec d > 2, la notation 5i(a, /3,7, (5) 
pour d = 1. 

La paramétrisation qu'on rencontre souvent dans les applications est celle 
de Samorodnitsky et Taqqu |20j . La différence entre notre paramétrisation 
et celle de |20| concerne la masse totale de la mesure spectrale a. Notons 
0"s le paramètre d'échelle si d = 1 et la mesure spectrale si d > 2 dans |20) . 
Dans le cas unidimensionnel 

7 = Cacrf, 

et dans le cas multidimensionnel 

Cj(-) = CaCTsi-), 



OÙ Ca est définie par (3.12). L'avantage de notre paramétrisation est que 



la mesure spectrale coïncide avec celle qui apparaît dans la propriété de 



régularité (2.4). Soit X la variable aléatoire de loi 5i(a, /3, 7, 5), alors on a 
r lim x"P{X > x} = o-(l), 

1 ""iST x°P{X < -x} =ct(-1), ^^'^^^ 

\ x—>-oo 

où 

(1 + /3) (1-/3) 

= — ^ — 7, o"(-l) = — n — 7- 



La démonstration de (3.13) peut être trouvée dans ^20j, page 16 et 197. 
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Propriété 1. Soit X le vecteur aléatoire de loi 5rf(a,(T, 5). Le vecteur aléa- 
toire X est StaS dans avec < a < 2 si et seulement si 

(a) a / 1, J = 0, 

(b) a = 1, Jgd-i s^^^a{ds) = 0, pour k = 1,2, . . . ,d. 

Propriété 2. Le vecteur aléatoire a-stable est symétrique si et seulement si 
ô = et a est une mesure symétrique sur S'^^^, i.e. 

a{A) = a{-A) (3.14) 

pour tout ensemble borélien A dans S"^^^. 



Données a-stables simulées 

Il existe des méthodes connues pour simuler les vecteurs aléatoires stables. 
Chambers et al. [2] ont proposé une méthode pour simuler les variables aléa- 
toires stables unidimensionnelles arbitraires. Plus récemment une technique 
précise a été construite dans |23| en utilisant une transformation non-linéaire 
d'une paire des variables aléatoires indépendantes de loi uniforme et de loi 
exponentielle. Une méthode pour simuler les vecteurs aléatoires stables basée 
sur la discrétisation de la mesure spectrale a été présentée par Modarres et 
Nolan |10) . Le résultat suivant présenté par R. LePage et al. montre qu'on 
peut simuler des lois stables en utilisant la représentation de série. 



Théorème 3.2. (I9l, Th. 3) La série définie par (1.2) avec a < 1 ou a 
symétrique a la loi StaS. 

Soit ^ Q, (j un vecteur aléatoire stable d'indice a et de mesure spectrale a 
vérifiant la condition du théorème précédent. On simule ^a,a en utilisant la 



somme partielle de la série de LePage (1.2) 



4 = J^cPri/^e,. (3.15) 

i=l 

Dans un premier temps deux lois unidimensionnelles 5i (0.75, 0.5, 1, 0) 
(non symétrique) et 5i(0.75, 0, 1,0) (symétrique) sont considérées. On simule 
les variables aléatoires en utilisant la somme ( |3.15 ) avec le nombre de termes 



k différent. Les valeurs 50, 500 et 5000 sont utilisées pour simuler la loi stable 
non symétrique, tandis que la loi stable symétrique est simulée avec k = 5,10 
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et 50. On observe que l'approche de la densité empirique à la densité réelle 
est plus rapide pour la loi stable symétrique, voir la figure [3] Cela signifie 
que la somme partielle (3.15) donne une bonne approximation pour les lois 



a-stables symétriques. Dans le cas non-symétrique, la convergence est lente 
pour une valeur de k modeste. 




— densité réelle 

■ ■ ■ k=50 

- . k=500 

■ - k=5000 



25 - densité réelle 

■ ■ ■ ■ k=5 

k=10 
k=50 



-5 -4 -3 -2 -1 1 2 3 4 5 

Figure 3 - Densités des lois 5i(0.75, 0.5, 1, 0) (haut) et 5i(0.75, 0, 1,0) (bas) 
simulées par la somme partielle de série de LePage avec différentes valeurs 
de k. Le nombre de données simulées est lO"^. 

Ensuite les vecteurs aléatoires 2-dimensionnels de lois 0.75 et 1.5-stables 
avec la mesure spectrale uniforme sont simulés. Le nombre de termes de la 
somme partielle est 10. On calcule la densité réelle en utilisant la fonction 
Matlab "mvstablepdf" dans le package appelé "STABLE" obtenu sur de- 
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mande au site http: / /www.RobustAnalysis.com/[ Les courbes de niveau des 
densités réelles et empiriques sont présentées dans la figure [4] 

Enfin, on simule les vecteurs aléatoires 2-dimensionnels de lois 1.5-stables 
avec les mesures spectrales ai, (J2 et dont les densités sont définies par 

•' ^ ' [0, smon, ^ ' 

..m = / ^œs(20)' ^G[-f,f)U[f,f), 

' \ 0, sinon, ^ ' 

fm = J|cos(20)|, 0G[O,2^). (3.18) 



Ces trois mesures vérifient la condition (3.14), les lois strictement stables 
correspondantes sont donc symétriques. Le package de programmes qu'on 
utilise ne permet pas de calculer la densité de la loi stable avec la mesure 
spectrale absolument continue. Nous présentons dans la figure [5] les densités 
de mesure spectrale fi et les courbes de niveau des densités empiriques de 
lois 52(1.5, (Tj, 0), i = 1, 2, 3. 

Un des usages de ces données est de tester la robustesse des procédures 
statistiques multivariées : engendrer aléatoirement l'ensemble de données et 
évaluer les statistiques par la procédure considérée. Un autre usage possible 
est de calculer la probabilité P{X E A} pour l'ensemble A <ZW^ on X est 
le vecteur aléatoire qu'on simule. En général le calcul numérique de cette 
probabilité est difficile, voir par exemple |14j et [15]. On peut les estimer 
par la méthode standard de Monte-Carlo une fois que nous savons générer 
les vecteurs aléatoires de loi prescrite. Le calcul devient coûteux si on prend 
le nombre de termes k grand. Donc cette méthode n'est pas avantageuse 
pour la simulation des lois stables non symétriques. Le troisième usage de 
ces vecteurs est la simulation dans les domaines d'application, par exemple 
l'analyse du portefeuille stable p^ . 



4 Conclusion 

Sous l'hypothèse de variation régulière tous les éléments StaSs avec a G 
(0, 1) peuvent être représentés par la somme des points d'un processus ponc- 
tuel poissonnien, dite représentation de LePage. La représentation de LePage 
de la loi strictement max-stable implique que la loi ne dépend que des points 
les plus grands du processus ponctuel poissonnien correspondant. En utili- 
sant cette propriété, on peut déduire les lois jointes et marginales des vecteurs 
aléatoires max-stables avec les mesures spectrales données. La propriété de 
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Figure 4 - Courbes de niveau des densités réelles (gauche) et empiriques 
(droite) de lois a-stables bivariées de mesure spectrale uniforme, a = 0.75 
(a.), 1.5 (b.), k = 10. Le nombre de données simulées est 10^. 
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variation régulière de la queue des lois est vérifiée. Les éléments aléatoires 
StaSs vérifiant la condition de variation régulière peuvent être simulés en 
utilisant la somme partielle de la série de LePage. Le temps de calcul de 
simulation est faible pour les lois strictement max-stables et StaSs symé- 
triques. Les résultats de simulation montrent visuellement que l'indice a et 
la mesure spectrale a mesurent respectivement la disposition des masses et 
la disposition directionnelle de probabilité autour du centre. 
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